i in h i 






COMPARISON OF METHODS 
IN CALCULATING FREQUENCIES 
OF CORNER-SUPPORTED 
RECTANGULAR PLATES 

by Robert E. Reed, Jr. 

Ames Research Center 
Moffett Field, Calif. 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. • SEPTEMBER 196S 



TECH LIBRARY KAFB, NM 

.... .am lllll HUB Hill Hill ilil III! \ i 'i-1 'Jlj , J 



TECH LIBRARY KAFB, NM 



0071113 

NASA TN D-3030 


COMPARISON OF METHODS IN CALCULATING FREQUENCIES OF 

CORNER- SUPPORTED RECTANGULAR PLATES 

By Robert E. Reed, Jr. 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 




COMPARISON OF METHODS IN CALCULATING FREQUENCIES OF 


CORNER -SUPPORTED RECTANGULAR PLATES 

By Robert E . Reed , Jr . 

Ames Research Center 


SUMMARY 


The lower natural frequencies of a rectangular plate that is pin supported 
at the corners are found analytically by two approximate methods . The first 
is the Ritz method and the second is a series solution to the differential 
equation of motion. These frequencies are also found experimentally for com- 
parison. All of the results are presented along with some previous work done 
by other authors using a finite difference method. Some of the advantages and 
disadvantages of the two methods pertaining to the numerical computations are 
discussed. 


INTRODUCTION 


The linear equation of motion, based on the Kirchoff plate theory, for 
the free transverse vibration of thin rectangular plates is well known but 
exact solutions have been found for just a few types of boundary conditions. 
For more difficult cases, approximate methods are available for calculating 
the natural frequencies and mode shapes and have been extensively applied by 
various authors to problems in which the plate is supported along entire 
boundaries (see, e.g., ref. 1, ch. 7)* In the following, two approximate 
methods, which can be applied to many other problems, are used for determining 
the natural frequencies and mode shapes of a rectangular plate with isolated 
pin supports at the corners. This problem is of interest since it is common 
design practice to support panels at isolated points (e.g., solar panels, 
hatches, covers, etc.) and also because the solutions point out some of the 
difficulties that can be encountered in plate vibration problems when 
approximate methods are used. 

The first method used is the well-known Ritz method. Natural frequencies 
and mode shapes of the lower modes were calculated to what appeared to be 
about 1 -percent accuracy by using more and more terms of the assumed series 
until the result indicated (at least to the author) the convergence to the 
desired ac curacy . However, some of the numerical examples calculated herein 
were previously treated by Cox and Boxer (ref. 2) who used the approximate 
method of solving the finite difference equations of motion on a digital com- 
puter. Upon comparison with their work (they felt their answers were accurate 
to within about 1 percent), discrepancies in natural frequency were found to 
be as large as 10 percent. Instead of exploring the convergence of the Ritz 
method, a series solution to the differential equation of motion was obtained 
and the calculated frequencies showed that, in most cases, the discrepancies 



were due to the slow convergence of the Ritz method. A short experimental 
program was then undertaken to see the correlation between experiment and 
theory. All of the results are presented herein for comparison along with a 
discussion of the problems associated with the methods and their adaptation to 
the computer . 


ROTATION 


Ani 

a 

®ni 

b 

D 

E 

g 

h 


•^o 

^max 

t 

^max 
W(x,y ,t ) 

w(x,y) 




% 

7 

A 


integration constants, i = 1, 2, 3, k 
length of plate 

integration constants, i = 1, 2, 3, k 
width of plate 

Eh 3 

bending rigidity of plate, — 

12(1 - V 2 ) 


Young f s modulus 
acceleration of gravity 
plate thickness 

summation indices 

a 

b 

amplitude of kinetic energy 
time 

amplitude of strain energy 
transverse displacement 
amplitude of W(x,y,t) 

Cartesian coordinates 

nrt 

a 

nit 

b 


weight density 

dimensionless frequency parameter, 



2 



V 


Poisson* s ratio 


a) circular frequency 

METHODS OF ANALYZING PLATES ON ISOLATED SUPPORTS 


The form of an approximate solution for a plate on isolated supports can 
sometimes be obtained from the solution for the more common problem of a plate 
supported along entire edges. For a rectangular plate., this latter solution 
can usually be chosen as 

W = m n a tnn F m( x ) G n(y) 

where the displacement W(x,y,t) = w(x,y) sin wt and the function F m (x)G n (y) 
and its derivatives can satisfy certain conditions along the lines 
x = constant and y = constant but cannot satisfy restraint conditions at 
isolated points. For a plate on isolated supports, a superposition of solu- 
tions of this type can be used so that each part satisfies conditions along 
certain, but different, lines and the intersection of these lines are at the 
isolated supports. For example, the sketch below shows a plate on three 
supports. Let w be of the form 



y Corner -Supported Rectangular Plates 

= ~ ~ 1 The plate being considered in the 

present report is shown in the adjacent 
sketch. The solution can be taken as 

b 

w(x,y) = vi(x,y) + w 2 (x,y) 

° where 


Corner -supported rectangular plate 


wi(0,y) = wa(a,y) = 0 
■w 2 (x,0) = w 2 (x,t>) = 0 
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The complete boundary conditions are: 


Displacements : 


Moments : 


w 


(0,0) = w(0,b) = w(a,0) = w(a,b) = 0 


^ (0,y) + v (0,y) = 


0 




dy 


S-l (a,y) + v Si (a,y) - 


dx 


Sy" 


w 


(XyO) + V (XyO) = 0 


dy^ c)x 

Effective shear forces: 


S-l U,b) + v Si (X,b) 
by Sx 2 


= 0 


b 3 w 

(0,y) 
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^ (0,y) = 0 

d 3 w 
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Sx 3 



bx Gy 
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The reactions at the corners are given by 2Mxy evaluated at the corners (this 
is the corner force that arises when the shear Q x or %■ is combined with the 
rate of change of the twisting moment M^y to give the effective shear V x 
or V y )(ref - 3) • 

Ritz method soluti on. - The Ritz method is well known (ref. *0 and is 
based on the fact that for a freely vibrating system, for which the motion is 
harmonic, the sum of potential and kinetic energy is a constant and the maxi- 
mum values (with respect to time) of these forms of energy are equal to each 
other. If the time dependent displacement is denoted by W(x,y,t), it can be 
written in the form 


W(x,y,t) = w(x,y)sin cot 


If w(x,y) is of the form 


w(x,y) 


a nrn^m l) (x)G^ l) (y) +b mn F( 2) (x)G( 2) (y) 


m n 


( 1 ) 
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the ratios of the unknown a 1 s and b ! s and the natural frequencies can be 
determined from the set of linear algebraic equations obtained from 




( v max " T max) - 0 
(^max “ T max) “ 0 
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where , in this case. 
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( 2 ) 


( 3 ) 


T max = ^ ^ < U > 

Each term of the series in equation (l) must satisfy the boundary conditions 
on the displacements and rotations but does not have to satisfy those on 
stress (ref. 4 ). Of course, one would expect that fewer terms of equation (l) 
would be needed to represent a mode if some or all of the stress boundary 
conditions were satisfied. Also, it is known that the frequency calculated 
for a particular mode from an approximate form of w will be above the exact 
answer (ref. k) . 

For the problem being considered, the displacement is chosen to be 
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+ V cos 
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( 5 ) 


where the single summation contains the m = 0 terms of the double summation. 
One should note that for each value of m and n, equation ( 5 ) satisfies the 
displacement boundary conditions but does not satisfy either the moment or 
transverse shear boundary condition. Equations (2) will be of the form 
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( 7 ) 


^max ^max 


Sa ij 

Saij 

^max 

^max 

db Q j 

8b 0 j 

^max 

^max 


= 0 


= 0 


_ © i — • • • 

db n - j 8b,- 4 0 = • • . 

1 J x J 


(8) 

( 9 ) 


Substitution of equation ( 5 ) and. its derivatives into equations (3) and 
( 4 ) and integrating over the surface of the plate will give V max and T max in 
terms of the a^'s and b-y's. Then, substitution of V max and T max into 

equations (6) through (9) will give the following algebraic equations from 
which the frequency parameter A 2 , where A 2 = (a 4 /* 4 ) (7hw 2 /Dg) and L Q = a/b, 
can be determined: 

a Qj2 (jV-A 2 ) + >j" l] I ^ [(-D n -1](^ 2 J 2 L 0 2 -^) 

n=i 


^ I ^ [ ( "l)^( ( m2 ^ 2 ^o 4 + vnS 3 2 ^ J 0 2 ~^ 2 ') ~ 0 (10) 


m,n=i 

b . 2 (j 4 -A 2 ) V ^22. [(-l) n -l](vn 2 j 2 L 0 2 -A 2 ) 

OJ Jt j i—i n 


^ I n(j5% a ) [(-D n -in(-l) J (-l) m -l](m 2 j 2 + vn 3 j a L 0 a - Ag ) = 0 (ll) 


m ,n=i 


ft J 


[ ( -l) n ( -l) 1 -l] (n 2 i 2 + vn 2 j 2 L 2 -A 2 ) 


n=i 

n^i 


^ V bmn n C(-l) 1 (-l) n -l][(-l) 3 (-D ,a -l] s±s +m a a 4 + (2 . v)a ^ , + m ^ L . 0 ( 12 ) 

* m,n=i (n 2 -i 2 ) (j 2 -m 2 ) 0 0 0 

m^j ,n^i 


2 ) 2 -* 2 ] + itLL g) 3 - 1 ] y [(-D^-D 1 -!]^ 2 ! 2 !, 
Si 

+ ly_ ^ a^[(-l) m (-l)J- 1 ][(- 1 )n(. 1 )i. 1 j 


to + vn 2 j 2 L 0 2 -A 2 ) 


m,n=i 
mfc) ,nj*i 


(n 2 -i 2 )(j 2 -m 2 ) 


[m 2 j 2 + n 2 i 2 L 0 4 + (2-v)m 2 i% 2 + vn 2 j 2 L 0 2 -A 2 ] =0 (l 3 ) 
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In equations (10) through (13), each value of i and j gives a single 
equation and m and n are dummy indices. The summation notation is 


I -E E<— 


m,n=i m=i n=i 
,n^i m^j n^i 


Inspection of these equations shows that they uncouple into four groups of 
equations. Each of these groups represents modes containing a certain type of 
symmetry or antisymmetry as follows: 


modes symmetric about the 


( 1 ) 
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= odd 
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= odd 
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a 03 * a 21^ * 

• • ^ 01 ; 
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• • 
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a. • : 

i = odd, j 
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^mn^ 

m = even ; n 

= even 





a ii' 

a i3 > a 31 ; * 
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• • 
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m = odd, n 

= odd 
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a 04 ^ a 24^ * 

• • b ii> 

^13* 
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m • 

w 

a ij > 

i = odd, j 
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a i4 2 a 3 2.* * 

• * ^ 12 > 
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lines x = — and y = §• 
2 J 2 


modes symmetric about y = ^ 
and antisymmetric about x = 


3/ 

modes symmetric about x = ^ 
and antisymmetric about y = 


modes antisymmetric about the 

a 
2 


lines x — and y = 


Combinations (2) and (3) are equal for a square plate and are related to each 
other by the ratio L 0 for rectangular plates. These four combinations of 
equations are given in appendix A. The equations of each combination can be 
treated separately and will have the matrix form 
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where the matrix is symmetric and all of the elements of the C array are 
nonzero. Each nonzero element of the matrix is of the form 

Although these equations are not in the standard eigenvalue form of 

[u] - = o 

the roots of the determinant of the coefficient matrix and the mode shapes can 
be found with existing subroutines . 

Trigonometric series so lution.- An appropriate solution to the differen- 
tial equation of motion can be found in the form of an infinite series. The 
differential equation is 


+ 2 - - £ 4 w = 0 
8x 4 8x 2 by 2 8y 4 

where 

e 4 = 

Dg 


(1^) 

( 15 ) 


Using the same method discussed earlier, one finds an appropriate form of 
w which satisfies the displacement boundary conditions as being 


oo 00 

W =X f n(y) sin + ^ Sn( x ) sin ^ (l6) 

n=i n=i 


Upon separation of variables, the unknown 
equations: 


a 4 f n (y) 

ay 4 


2 a 2 f n (y) 
2a n — V- + 
ay 


functions must satisfy the following 
(c^ 4 - £ 4 )f n (y) = 0 (17) 




a 2 g n (x) 

dx 2 


+ 



^ 4 )g n (x) = 0 


(18) 


where oc^ = mt/a, P n = mt/b; and f n (y)sin(n«x/a) is assumed to be linearly 
independent of g n (x)sin(nrty/b) . Some confusion can result if they are not 
independent, that is, f n (y) = A sin(njty/b) and g n (x) = B sin(mtx/a) . 

Appendix B discusses this problem. 

Since the coefficients of equations (17) and (18) are constant, the solutions 
are easily found to be 

f n (y) = A ni cosh 0 n y + A n2 sinh 0 n y + A n3 cos 0 n y + A n4 sin 0 n y (19) 

g n (x) = B ni cosh cp n x + Bj^g sinh cp n x + B n3 cos cp n x + B^^. sin cp n x (20) 
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•where 


0 n = + «n 2 0 n = A 2 - “n 2 

<P n = + p n 2 *Pn = ~ P n 2 








( 21 ) 


It should he pointed out that can rea ^- or imaginary; depending on 

the relative value of A and ot^ or (B n . If % or cp n is imaginary, then the 

trigonometric function in equations ( 19 ) or (20) becomes the corresponding 
hyperbolic function. 


Referring to the boundary conditions given by equation ( 3 ), equation (l6) 
is seen to identically satisfy the zero displacement conditions at the corners. 
The conditions for zero bending moment give the following: 


a 2 g n (°) 


ax" 


a g n (a) 


dx 


d 2 f n (0) 

<iy 2 

a 2 f n (b) 

ay 2 


v Pn 8n(°) = 0 
2 . 

v(B n g n (a) = 0 
V(x^ fji(O) = 0 
v % 2f >> 0 0 




( 22 ) 


J 


The condition of the effective shear being zero gives the following equations: 


i -» [t-> ^ - %%«] * 1 - <->%■ ^ 


( 0 ) 




sin P n y = 0 


n=i 


00 00 

X aj-l ) n [(2-v) ^2 - «n%(y)] + X - ( 2 ' v)|3 n 2 sin ^ = 0 

n=x n=i 

OO CO 

V D r, , d 2 g n (x) 2 1 y rd 3 f n ( 0 ) df n (o)l . 

1 Pn [(2-v) - 5 ^- - P n SnWj + ^ ( Hy-] s “ % x = 0 

n=i n=i 

V , ,n T, x d \( x > 2 . ,1 f 1 ra 3 f n (b) , x 2 

1 P n (- d ) [(2-v) - P n g n (x)J + ^ [ - (2-v)a n sin a^x = 0 


> ( 23 ) 


J 
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For the bending moment, equations (22) show that the constants can he deter- 
mined so that each term of the series in equation (l6) identica lly satisfies 
the condition of zero moment on the edge. Equations (23) show that the effec- 
tive shear boundary condition is satisfied only by a sum of terms. Equa- 
tions (23) can be simplified by multiplying by the appropriate orthogonal 
trigonometric function and by integrating over the plate interval. For 
example, the first one of equations (23) becomes 



a n 2 f n (y) 


sin f^y dy 


+ b l~ d3 gm(°) 
2 L dx 3 


(a-Wiy 2 


agJ 0 )] 

dx 


= 0 


(210 

The integrals in equation (24) can be evaluated and the resulting four sets of 
equations are in terms of the eight sets of -unknown constants, A ni , . . ., An4> 
Bn!, . . B n4 . For the integrated equations it is assumed that / Pm* 

The condition that 0 n ^ is connected with the assumption that 

f n (y)sin o&qx is independent of g m (x)sin l^y and is discussed in appendix B. 


Four of the eight sets of constants can easily be eliminated from 
equation (24) by use of equations (22). This results in having a set of 4 n 
simultaneous homogeneous algebraic equations to solve (after truncating the 
infinite series at n = N terms) . These equations are given here because 
they can easily be programmed for a digital computer, and additional modes and 
frequencies, to those given herein, can be calculated. (See the Numerical 
Computations section for a discussion of some of the numerical problems and 
ways to, at least partially, avoid them.) In these equations, the following 
notation is used: 

A 2 = f 4 sl - a4yhc ° 2 
it 4 rt 4 Dg 


0 n* - £ A - " a > 
cp n * = Jt - n 2 L Q 2 * 


= £ 


+ n c 


9. 


> n * = n n/a + n 2 L Q 2 


(25) 


2[(2-v)A 3 + ( l-v) 2 | 2 p 2 A] 

[ -A 2 + U s +p 2 ) 2 ][A + (l-v)i 2 ] 


( 26 ) 


P (|) [A - (1 -vH 2 ] 2 

[A + (l-v)i 2 ] 


(27) 
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P 3 U) ='/7'-i 2 [A + (1-V)5 2 ) (28) 

where £,p are used here as dummy variables. The equations are 
N 

^ n(-l) n Pi(n,mL 0 ){-A tl3 [l - (-l) m cos 0*] + (-l) m A n4 sin 0 n *} 
n=i 


it L -n t r \ 00 8 cosh cpm* - 1 „ , T x 

Pa^o) — + P s (mL 0 ) 

L L sinh cftn* 


sin cp 


* 


m 


Pm4 


P 2 (mL 0 ) 


cosh cp m * sin cp m * 


sinh c^* 


P 3 (m L o) c os cp m * j- = ° ( 29 ) 


N 


nP 1 (n^mL 0 ){-A n3 [l - (-l) m cos 0 n *] + sin 0 n *} 


n=i 


2 mL^ 


+ 0 T 2 I ^3^2 


cos - cosh cpm* 

sinh cp * 
m 


+ B, 


m4 


P 2 (mL 0 )sin cp * 


sinh cp * 


m -p 3 (^ 0 ) 


= 0 (30) 


2mL- 


2 }^ m 3 


P 2 (m) 


cos 0 m * cosh — 1 

sinh 0 * 
m 


+ Ps(m)sin 6^* 


+ ^14 


P 2 (m) 
N 


cosh 0 m * sin 0 m * 
sinh 0 m * 


P 3 (m)cos 0 m * 


^ n(-l) n P 1 (nL 0 ,m) - (-l) m cos cp n *] + (-l)" 1 ^ sin cp n *| =0 (31) 

n=i 
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cos e m * - cosh e n * 
sinh 0 m * 


+ ^m^ 


P 2 (m) 


sin 0 m * 


sinh 9, 


iv - p s w] } 

m J J 


N 

^ nP 1 (nL 0 ^m) - (-l) m cos cpn*] + (-l) m B n4 sin 9 n *j" = 0 (32) 


n=i 


where A f n 2 + m 2 ]^ 2 or A f n 2 ]^ 2 + m 2 (see appendix B) . Each value of m 
gives a separate equation so the resulting coefficient matrix will be of order 
4N. The frequency parameter A is obtained from the condition that the 
determinant of the coefficient matrix is zero. One significant difference 
between these equations and those of the Ritz method is that there is no form 
of 


[A] - A 2 [B] = 0 

to these equations since A is contained in 0 m *, %i*-> This difference 

is discussed in the next section. 


NUMERICAL COMPUTATIONS METHOD 


A discussion of the computation is given not because any new numerical 
scheme is used (the programs that have been written are made up of existing 
subroutines) but because each method of solution has certain numerical advan- 
tages and disadvantages. For the problem being considered, the Ritz method 
offers comparatively easy computations with possible low accuracy and, con- 
versely, the series solution offers high accuracy with possible serious numer- 
ical difficulties. So, depending on one's needs, either method may be the 
best to use for a given problem. 


Ritz Solution 

The form of the Ritz method equation is fairly convenient for computation 
on a digital computer. For the modes calculated, the matrices were small 
enough so that no major problems were encountered. The roots were found by 
the computer varying the frequency over some specified range and detecting 
when the determinant changed sign. An iteration scheme was then used to con- 
verge on the root. The computing was done in single precision. Eventually, 
if larger matrices are used, a point would be reached where the accuracy had 
diminished to an unacceptable level but double precision could be used and 
one might, at least partially, normalize the equations. For instance, the 
constants a 0 j , b Q j , a^ , bpj could be replaced by a^/ j 2 L Q 2 , b^j/j 2 ^ 

aj_j/(i 2 + j 2 L Q 2 ), b^j/(i 2 L Q 2 + j 2 ), and the denominators retained in the 
coefficient matrix. 
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Series Solution 


The equations (29) through (32), as given, can present some serious 
computational difficulties. One problem is that A cannot be factored out in 
the form [A]-A 2 [b] = 0, so the standard eigenvalue routines cannot be used. 

The only feasible way of obtaining the frequencies seems to be to plot the 
determinant as a function of A and to note the zero values of the determi- 
nant. However, as the determinant approaches zero, accuracy is lost and there 
will be a range of frequency within which the computed determinant has zero- 
place accuracy and therefore is meaningless. If this range is larger than the 
range needed to interpolate the frequency to some desired accuracy, the fre- 
quency cannot be determined. In the present case, this problem became more 
prevalent as the plate’s length -to -width ratio or matrix size increased. In 
order to include all of the frequencies listed herein, it was necessary to 
rearrange equations (29) through (32). Comparing equations (29) and (30) and 
equations (31) and (32), one sees that the terms inside the summation are 
identical except for the term (-l) n . By adding or subtracting equations (29) 
and (30) for each value of m (the same applies to eqs. (31) and (32)), either 
the odd A ns , A n4 T s (n = 1,3*5>- • ■) or the even A n3 , A n4 ' s will vanish. 

This introduces many more zeros into the matrix and eases the loss -of -accuracy 
problem. These equations were then normalized by substituting the following 
for the given constants: 

An3 ^T14 

cosh ^ sinh — 

f'o Lo 

One should be careful to normalize with positive nonzero quantities since the 
quantity, if ever zero, will introduce a point of discontinuity, and possibly 
a change in sign, into the determinant. When the determinant is plotted, 
these points may appear to be natural frequencies. The reason for choosing 
the particular normalizing factors shown is that for large values of m 2 , say, 
m 2 > 30, and small values of A, say, A < 10, the large terms in the matrix 
are of the form 


^m3 ^K14 

cosh mot sinh m at 


cosh cp m * = cosh at \l~K-\m 2 = 0(cosh mat) 

cos 0 n * = cosi at \ 1 m 2 -A = cosh at \l m 2 -A = 0(cosh mat) 

so that normalizing these terms by the factor cosh mat makes them about 
unity. But if higher frequencies are desired such that A is not small com- 
pared to m, this factor may not be adequate and it should be weighted more 
heavily toward the effect of A. So, depending on the range of the parameters 
(m, L 0 , A, etc.) being considered, considerable trial-and -error rearrangement 
of the equations may be necessary. 

Finally, double precision was used in the computations and in generating 
the elements of the matrix. The revised equations are given in appendix C. 
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EXPERIMENTAL METHODS 


A short experimental program was undertaken to see the correlation with 
the theoretical results. Preliminary testing of a plate that was visibly 
warped indicated that the frequency was quite sensitive to plate warpage . 
Therefore, some care was taken to insure that the plates for which the results 
are presented were flat (a "flat” plate is one in which the maximum deviation 
is a small percentage of the thickness) . The method of support is shown in 
figures 1 and 2. The supports were adjusted until there was no loss of con- 
tact between the support and the plate corner during vibration. This means 
that some compression was introduced into the plate, but its effect was con- 
sidered negligible since the force was close to a point load which produced a 
very small average compressive stress throughout the plate. Changes in the 
support adjustment gave no measurable difference in natural frequency. 

Several methods of excitation were used. The primary one for measuring the 
frequencies was a 30 -watt sound speaker that excited the plate through an 
attached string (see fig. l) . The speaker was driven by a variable frequency 
oscillator and the motion of the plate was measured by a capacitance -type dis- 
placement instrument. Since the string transmitted no net compression, the 
exciting force was not a complete sine curve but this was considered to have a 
negligible effect on the measured natural frequency. 

Figures 3 and k are examples of the often used method of detecting the 
mode shapes. Sand (carborundum chips were used for the pictures because of 
their dark color) is sprinkled on the plate during vibration and the particles 
migrate to the node lines. 


RESULTS AND COMPARISON 


In the following remarks, the modes are referred to by number (see 
table I) . This can be misleading since the relative magnitude of the frequen- 
cies associated with two mode shapes can change as L 0 changes (see modes 6 
and 7, table I) . But for reference purposes, the modes are numbered in the 
order in which the frequency occurs over the major portion of the range of 
L q considered. 

Table I is the summary of the calculated frequencies along with the 
corresponding ones of reference 2. The results obtained from the Ritz method 
were calculated from the mode shapes given in table III, where the largest 
matrix considered was 12 X 12. The series solution used matrices of order 1 6 , 
20, and 2k and all three different sized matrices gave the listed value of 
frequency for most of the modes and a maximum variation of two digits in the 
last significant figure for the other modes. The series solution should be 
considerably more accurate for a given matrix size than the Ritz method 
because each term of the series satisfies the differential equation of motion, 
the displacement and moment boundary conditions. The comparison between this 
solution and reference 2 is seen to be very good. The frequencies in refer- 
ence 2 were obtained by extrapolating the frequencies found by considering the 
plate divided into a square gridwork in which the width of the nonsquare 



plates ¥as divided into k and 6 grid spacings and the square plate* s width was 
divided into 5 and 6 spacings. Richardson 1 s h 2 extrapolation formula 
(ref. 2) was used to obtain the value for an infinite number of mesh points. 

There are several interesting points that can be seen from table I. 

Mode 1 approaches the fundamental frequency of a simply supported beam as Lq 
increases. The beam frequency is A = 0.95^ (the beam frequency is adjusted 
since A, as defined here,, contains Poisson's ratio) , and A = O. 95 I for the 
plate when L 0 = 2.5* Modes 2 and 3 for Lq = 1*0 have the same frequency but 
different mode shapes (although, considering symmetry, they are the same) . 

This occurrence is not uncommon in plate problems and at that frequency the 
mode shape can be either one or some superposition of the two mode shapes. 

Mode k- is the same as a flexural mode of an unsupported plate . Ritz treated 
the case for a square plate and the corresponding frequency parameter is 
A = 2.00 (Ritz's calculations, reported in ref. were for v = 0.225^ so the 
value was scaled up to v = 0.3 by use of ref. 2). It is interesting to note 
how this mode shape changes for a nonsquare plate (see table I and fig. 3 ) • 

In figure 3, it appears as if the node line does not intersect the support 
point but theoretical results show that is passes very close to the edge of 
the plate in the region of the corner. Later tests showed this discrepancy 
was caused by some slight support motion due to the vibration of the test fix- 
ture which apparently had a natural frequency near the exciting frequency. 

It is also of interest to see the number of other problems solved. 

Looking at the mode shapes in table I, all of the straight node lines are 
lines of antisymmetry along which the bending moment is zero. These, then, 
correspond to simply supported boundaries and can be considered to be plate 
edges. So, for example, modes 2, 5; and 6 give the first, second, and third 
mode of a plate simply supported along one edge and point supported at the 
other two corners (for the range of L 0 given and substituting, of course, 
the proper length and width of the plate) . Mode 4, for L 0 = 1.0, gives the 
fundamental mode of a ^5° right triangular plate simply supported along two 
edges, or the second mode of a ^5° right triangular plate supported along the 
hypotenuse and at the opposite corner . 

After the results were compared, the Ritz method was extended for modes 
1, k, and 7 for Lq = 1.0 to see how slow the convergence was. This result is 
shown in figure 5 . Mode k is seen to converge very fast. For just two terms, 
which are 

w = a 01 (^sin - sin 

the answer is within 1 percent. However, mode 7 has an error of about k per- 
cent after 32 terms. For modes 1 and 7> the results for matrices of order 2, 
8, 10, and 12 suggest that the series has converged (see dashed lines in 
fig. 5) when, actually, convergence has not been reached even with a matrix 
of order 32. This shows that even though one knows an approximate frequency 
is an upperbound, it is hard to tell the accuracy of an answer without 
comparison with a known solution. 



15 



Table II shows the comparison between experiment and theory (the series 
solution and ref. 2 ) and the correlation is very good. No attempt will be 
made to explain the small differences because the experiment was of a careful, 
but not precise, nature. 


CONCLUSIONS 


One important conclusion common to many plate problems is that one has to 
be very careful in tracing the modes as the plate changes its dimensions (i.e., 
from a square to a rectangle) . The relative order of frequency can change as 
can the mode shapes. The fact that mode shapes can change is especially 
important when one is concerned about the location of nodal and antinodal 
regions . 

The Ritz method yields equations which are easily adapted to the computer 
but the convergence is hard to predict . The fact that the frequency is about 
the same for different numbers of terms hopefully means that the assumed dis- 
placement has converged to the correct answer, but it may mean that the 
assumed displacement is converging very slowly and many more terms are needed. 
Of course, one would expect a higher accuracy for a given number of terms if a 
displacement function that satisfied at least some of the stress boundary 
conditions could have been chosen. 

The series solution presents some serious computational difficulties for 
certain ranges of the parameters but for this particular problem one has more 
confidence in its accuracy for a given matrix size since all of the boundary 
conditions are satisfied except the effective shear condition by each term in 
the series. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., June 28 , 1965 
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APPENDIX A 


UNCOUPLED EQUATIONS OF THE RITZ METHOD 


It was pointed out in the text that the sets of equations for determining 
the frequencies by the Ritz method uncouple into four groups of equations. 

The equations of each group represent modes with certain symmetry character- 
istics ^ so if a certain class of modes is desired, considerable simplification 
is obtained if the appropriate group of equations is used. Using the notation: 


R^i^m) = n 2 L 0 2 (i 2 L 0 2 + vj 2 ) 

R 2 (i,j,n) = n 2 (i 2 + vj 2 L Q 2 ) 

R 3 (j,m,n) = j 2 Lo 2 (m 2 L D 2 + vn 2 ) 

R 4 (,j,m,n) = j 2 (m 2 + vn 2 L 0 2 ) 

R 5 (i,j,m,n) = n 2 i 2 + m 2 j 2 L Q 4 + (2-v)m 2 i 2 L 0 2 + vn 2 j 2 L Q 2 
R s (i,j,m,n) = n 2 i 2 L Q 4 + m 2 j 2 + (2-v)m 2 i 2 L 0 2 + vn 2 j 2 L 0 2 


the equations are: 
Combination (l) 

a o,j2( j 4 L Q 4 -A 2 ) + 1| 




I 


n=x,3 




b on^ v n 2 j 2L o 2 -A 2 ) 


n 


Jt2 


I 




m=2^4^ . . . 

• ■ • 


n(j 2 -m 2 ) 


= 0 


(Al) 


t> oJ 2(j4_A2) + 


i6_ 

rt 2 j 


I 


a on ( vn 2 j 2L o 2 -A 2 ) 


n=x,3 




n 


i§3. V 

Li 


n(j 2 -m 2 ) 


= 0 


(A2) 



|| | II II || || II I mi | | nil IIH 


a ij [ (i 2 +«j 2l t> 2 ) 2- ^ 2 l + 


16 

jt 2 j 


b on n[R 2 (i, 0 >n)-A ] 


n=i*3. 


+ 


2 .2 

n -i 


l6j 


z 


m=2 , 4 y . . . 
31=1 J 3 y • • • 


~ b mn n [ R 5( 1 ^j^m,n)-A 2 ] 

(n 2 -i 2 )(j 2 -m 2 ) 


= 0 (A3) 


VaV+R) 2 -^ 


n=i,3j 


a on n[Ri(i^j^m)-A 2 ] 

2 • 2 
n -i 


+ 


l6j ^ a^ntReCi^j^n) -A 2 ] 




m = 2 _ ? 4 < , . . . 

n= • • • 


(n 2 -i 2 )(j 2 -m 2 ) 


= 0 (A4) 


Combination (2) 


4 2 l6j r a^fluC^m^-A ] 

l3 oj 2 (j- A ) + — Z_, 


tT *— < n(j 2 -m 2 ) 

• • • 

H = 1^3 ^ • • • 


= o (A5) 


a 1J [(i 2 +3 2 L 0 2 ) 2 -X 2 ] + i|r 


z 


b on n[R 2 (i.,j,n)-A 2 ] 


2 . 2 

*"J n - 1 


l6j y ~ b m n n| ^ B 5( 1 ^^ I1 )- A2] 

rt 2 Zj (n 2 -i 2 )(j 2 -m 2 ) 

m = 2 j 4 y # • • 

n = 2 , 4 , . . . 


= o (a6) 


b ij[( 1 % 2+ «3 2 ) 2 -A 2 ^ X 


a^nlRsCijjjm^-A 2 ] = q (a?) 


in — l j 3 ^ • 

n = i , 3 , . . . 


Combination (3) 


a 0;3 2(j 4 L 0 4 -A 2 ) + 


l6j 


I 


(n 2 -i 2 )(j 2 -m 2 ) 


b mn[R3( j ^ m ^n)-A 2 ] 


^ i! — 1 n(j 2 -m 2 ) 

m = i ^ 3 ; ... 

y 3 y m • • 


Hj [ ( i2+ j 2L o 2 ) 2 - A2 3 + l^ Yj 


m = i , 3 , . . . 
n = i , 3 , . . . 


(n 2 -i 2 )( j 2 -m 2 ) 


= 0 (A8) 


= 0 (A9) 
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^ij[(i 2 L 0 2 +j 2 ) 2 -A 2 ] +4r V 

L — i 


2 

a n[Ri(i.,j,m)-A ] 


H — 2 J 4 y • • • 


n 2 -i 2 


y 

7^ <3 


a^ntRsCi^j ,m,n)-A 2 ] 


= 2 ,T, . . . (n 2 -i 2 )(j 2 -m 2 ) 


nr 

n=2 j 4 y • • 


Combination (^) 


.a ,a T 2,2 -, 2l , l6j V T3r nn n[R 5 ( i ; ,j,m,n)-A^] 

+J L ° > - A 1 + — L - Wk«- 


m= i j 3 j . • . 
n = 2^4 <? • » • 


•b i j[(i 2 L 0 2 +j 2 ) 2 -A 2 ] + 


i6j y 

Tt 2 L 


hi l j 3 j • • • 
n 2 y 4 ^ • . . 


a ^ntRsCi^j ^iu,n)-A 2 ] 
(n 2 -i 2 ) ( j 2 -m 2 ) 


0 (A10) 


0 (All) 


0 (A12) 
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APPENDIX B 


DISCUSSION OF ASSUMPTION IN THE DERIVATION OF EQUATION (l4) 

In the derivation of the solution to equation (l4) , it was assumed that 
f n (y)sin a^x was independent of g n (x)sin P n y and later, in integrating the 
equations for the shear boundary condition, it was assumed that © n ■/= m which 
is the same as assuming A ^ n 2 + m 2 Lo 2 . When the frequencies are found by 
plotting the determinant versus A , the points A = n 2 + m 2 L 0 2 cannot be 
calculated from the equations in their given form since some terms become 
indeterminant. But, assuming the determinant is a continuous function of A, 
it appears as if it is identically zero for all the matrix sizes considered 
for A = n 2 + m 2 Lo 2 . One recognizes that these frequencies are identical to 
those of a plate that is simply supported on all edges. If 


f n (y)sin ct^x = g m (x)sin P m y = a^ sin P m y sin c^x 
the displacement is 


w = a^ sin P m y sin c^x 
and is a solution to the equation of motion if 

A = n 2 + m 2 LQ 2 

But, in order to satisfy the shear boundary conditions, the following type of 
equations must hold: 


W^nK 2 " (2-v)P m 2 ]sin P m y = 0 
^lyPm 2 - (2-v)a n 2 ]sin ax = 0 

These equations are for zero shear at x = 0 and y = 0. It is seen that 
a 7nJ1 = 0 . Therefore , the plate cannot have the same frequency and mode shape 
as a simply supported plate, that is, it cannot have zero displacement and 
zero shear along the same edge without having zero displacement everywhere . 

The remaining question is whether the plate can have the same frequency but a 
different mode shape . An example of this is that a beam with built-in ends 
has the same frequencies but different mode shapes as a beam with free ends. 
This question does not appear to have an obvious answer, in general, but the 
equations were modified to allow A to have a value of, say, n 2 + m 2 L Q 2 , and 
it was apparent that considerable uncoupling occurs . (For instance , 
sin 0—* = sin cp_* = 0; 1 - (-l) m cos 0g-* = 0.) It is easily seen in the small 

matrices how the constants A n3 , A n4 , Bng, B^ all have to be zero. A 4 X 4 
and an 8 X 8 matrix were investigated and the equations uncoupled in such a 
way that all constants had to be zero so there was no motion at that 
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frequency. These frequencies are not present either in the Ritz solution or 
in the experimental work so it is concluded that even though they represent 
roots of the determinant , they do not represent a mode that has motion. 
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APPENDIX C 


MODIFIED EQUATIONS FOR THE SERIES SOLUTIONS 


The modified equations which avoid numerical problems encountered for 
some of the modes are given by 


±2nPi(n,mL 0 ) 


-Aiatl - (-l) m cos 0 n *] (-l) m Aj 14 sin 


n=2,4, . . 
11=1,3, . . 


cosh ^ 


. , ?tn 
sinh 7— 
ho 


P 2 (niLo) (cosh cpm* ± l)(cos c^* + l) P 3 (mL Q )sin 


sinh cosh mrt 


cosh nut 


~ P 2 (mL 0 )sin ^(cosh cftn* ± l) _ P 3 (mLo)(cos cfa* ± l) l| = Q j 

sinh cp * sinh nut sinh mot _!/ 


P 2 (m)(cosh 0 m * ± l)(cos 0 m * +l) P 3 (m)sin 0 m * 


sinh 0 m * cosh ^ 


, P 2 (m)sin 0 m *(cosh 0 m * ± l) P 3 (m)(cos 0 m * ± l) 
^ m4 . sinh 0 m * sinh ^ sinh ^ 


±2nP 1 (nL 0 ,m) 


•B' [l " ( _1 ) m cos c Pn*-l ( _1 ) mB n4 sin ^n 


cosh mt 


sinh n?t 


0 ( 03 , 4 ) 


n=2,4, . . . 

n= 1 ,3 , • • • 


where for 


equations (Cl,3) n = 2,4,6,, . . m = 1,2, 3*4,. . . 

± use +, T use - 

equations (C2,4) n = 1,3,5,* • • m = 1,2, 3,- 4,* • • 


4- use + 
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TABLE II.- COMPARISON OF 





ro 


VJl 


AND THEORI 


Plate 2 10 X 20 X 0.173 inch 

2024 aluminum 

E = 10 .6X10 6 , psi (hook value) 

I j 

Ref. 2 !■ Experiment Series solution 

38.8 cps 38.3 cps 38.8 cps 


113 115 
136 137 
214 218 
261 267 


294 ' 298 




ro 


TABLE III.- RITZ METHOD MODE SHAPES 


Mode 1 


Mode 2 



a/b 

Frequency 
parameter > 
A 

Normalized mode -shape 
coefficients 

Frequency 
parameter , 
A 


Normalized mode- 
coeff icients 

shape 



a oi 

a 

1.000 

^01 

a 

1.000 


a 02 

a 

- 0.1248 

*bn 

a 

1.000 



a 03 

a 

- .0663 

603 

= 

- .0663 


a 04 

a 

- .0075 

6 13 

a 

-.0671 

1.0 

0.756 

a 21 

= 

.1737 


a 

.1737 

1.702 

a 22 

a 

.1695 

^31 

a 

- .0574 



a 23 

a 

.0329 

^23 

a 

.0329 


a 24 

a 

- .0055 

t >33 

a 

.0348 



a 4 i 

a 

- .0267 

^41 

a 

- .0267 


a 42 

a 

- .0146 

651 

a 

- .OO83 



a 01 

a 

.0869 

boi 



1.000 


a 02 

a 

-.1753 

6n 

a 

1.000 



a 03 

a 

- .0150 

^03 

- 

- .0320 


a 0 4 

a 

.0002 

^13 

a 

-.0499 

1.5 

.933 

a 21 

a 

.0950 

"^21 

a 

.0281 

2.308 

a 22 

a 

.1530 

^31 

a 

- .0748 



a 23 

a 

.00 56 

^23 

a 

.0161 


a 2 4 

a 

-.0059 

^33 

a 

.0205 










a 42 

a 

.0012 

6 5 1 

a 

-.0062 

1 


a 41 

a 

- .0102 

641 

a 

- .0068 


a 06 

a 

.0009 

6 15 

a 

-.0065 ■ 



a oi 

- 

- .0054 

’box 

a 

1.000 


a 02 

a 

-.1915 

bn 

a 

1.000 



a 03 

= 

- .0052 

^03 

a 

- .0179 


a 04 

a 

.0067 

^13 

a 

-.o 46 o 1 










a 22 

a 

.1537 

■ b 31 

a 

-.0893 ! 

2.0 

•958 

a 21 

a 

• 570 


a 

.oo 46 

2.941 

a 24 

a 

- .0088 

^33 

a 

.0159 


1 ] 

a 23 

r = 

•0013 

t >23 

a 

.0080 


a 42 

a 

.0080 

b 5 x 

a 

- .0028 


! 

a 41 

a 

- .0035 

t> 4 l 

a 

- .0025 


a os 

a 

.0007 

b is 

a 

- .0093 


1 

a 01 

a 

-.0197 

■box 

- 

1.000 


a 02 

a 

-.1984 

611 

a 

1.000 



a 03 

a 

- .0020 

b >03 

a 

- .0108 


a 04 

a 

.0100 

bi 3 

a 

- .0451 










a 22 

= 

•1550 

t >31 

a 

-.0967 

2.5 

.961 

a 21 

a 

.0372 

b 2 i 

a 

- .ooi 4 

3 -6 o 

a 24 

a 

- .0103 

1>33 

a 

.0135 



a 23 

= 

.0002 

t >23 

a 



a 42 

a 

.0120 

651 

a 

- .0009 



a 41 

a 

- .0007 

1)41 

a 



a 06 

a 

.0005 

"bis 

a 

- .0113 







TABLE III.- RITZ METHOD 



1 Frequency 
a ,/b ' parameter , 
A 


Normalized mode -shape 
coefficients 


a xl = 1.000 ; b 0 2 = - 0.1248 

a i3 = -.0671 b 04 = -.0075 

1.0 1.702 a 3 i = -.0574 b 22 = -1695 

a 33 = .0348 b 2 4 = — .0095 

a 51 = -.0083 b 42 = -.0146 

' * 

an = 1.000 ' b Q2 = -1539 

a i 3 = -.0817 b 04 = -.0380 

i 1.5 2.811 a 3i = .0205 b 22 = .1850 

; a 33 = .0428 b 2 4 = .0115 

1 a 51 = -.0209 b 42 = --0230 

an = .7924 b 02 = 1.000 

a i3 = -.0877 b 04 = -.0808 

2.0 3-52 asi = .1496 b 22 = .1713 

a 33 = .0426 b 2 4 = .0334 

a 51 = -.0304 b 42 = -.0272 

an = .1928 b Q2 = 1.000 

a i 3 = -.0337 b 04 = --0479 

2.5 3-77 asi = .1130 b 22 = .0536 

a 33 = .0150 b 2 4 - .0191 

a 5i = --0109 b 42 = -.0112 


ro 
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SHAPES - Continued 


Mode 4 



Frequency 
parameter , 
A 


Normalized mode -shape 
coefficients 


a Q1 = 1.000 ' b Q1 = -1.000 

aQ^ = —.0244 b 03 = .0244 

1.986 ^ a 2 i = -.0802 b 2i = .0802 

a 23 = .0112 b 23 = - .0112 


a 4 i = .0049 

b 4 i = 

-.0049 

1 

a 0 i = 1.000 

^01 = 

-.8108 1 

a 03 = -.0536 

^03 = 

.0002 

a 2i = .0756 

^>21 - 

.1693 

a 23 = -0255 

bs3 = 

.0171 

a4i = - .0107 

b41 = 

-.0191 

a 01 = 1.000 

boi = 

-.742 


a 03 = -*0726 b 03 = .0903 

5.69 a 2 i = .2048 b 2i = .2277 

a 23 = .0322 b 23 = .0341 

■ a 4 i = .0011 b 4 i = - .0277 

a 0i = .8438 b 01 = -.8072 

a.03 = -.3091 b 03 = .1068 

10.30 a 2i = -.8536 b 2 i = 1.000 

a 23 = .2419 b 23 — -.1570 

a 4 i = .0530 b 4 i = -.1440 




TABLE III.- RITZ METHOD MODE SHAPES - Continued 


-( 


a/b 

Frequency 
parameter , 
A 

Normalized mode -shape 
coefficients 

Frequency 
parameter , 
A 

Normalized mode -shape 
coefficients 

1.0 

4.20 

Q-12 *“ 1*000 
a 14 = -.1458 
0-32 “ *2107 

a34 = .06^-5 

a 52 = -.0433 

a 16 = --0O53 

bi 2 — 1*000 
bi4 = -.1458 
b 32 = -2107 

b34 = .0845 

T3 52 = -- 0 ^33 
bi 6 = -.0053 

5-23 

a n - 1*000 
a 13 “ “.0629 
a 3l = -.2789 
a 33 “ -.0007 

a 51 = .0318 

1=02 = -0.9798 
b 04 = .0944 

^22 = .0337 

b24 “ ” .0492 

b42 — - .0000 

1-5 

5.67 

ai 2 = - .0666 

a 14 = -.0342 

a 32 = • 1809 

a 34 = .0119 

a 52 = --0049 
a i6 = .0014 

bi2 = 1.000 

bi4 = -.0797 
b 32 = --0153 

b34 = .0328 

b 52 = -.0190 

bis = -.0101 

5-85 

an = 1-000 

a 13 = -.0446 
a 3 i = -.1487 

a 3 3 = .0225 

a^j_ = .0084 

b 02 = -.8072 
b 04 = .0477 

b 22 = .1222 

b 2 4 = -.0232 

b 42 = -.0078 

2.0 

6 .80 

a 12 = - • 1869 

a 14 = -.0122 

a3 2 = • 1784 

a34 = .0011 

a 52 = .0088 

a i6 = .0015 ’ 

bi2 = 1*000 
bi4 = -.0671 
b 32 = -.0548 
b 34 = .0238 

b 52 = -.0114 
bis = -.0134 

7-40 

a n = 1*000 
a 13 = -*0597 

a3i = -.0626 
^33 = .0316 
a 51 = “*0013 

b 02 = -.8479 
b 04 = .0324 

t>22 = .1557 

b 2 4 = - .0051 
b 42 = -.0151 

2-5 

7.95 

1 

a 12 = -.2238 
ai4 — - .0021 

a 32 = *1772 

^34 = -.0038 

a 52 = *013^- 

a l6 = *0013 ! 

b^2 = 1*000 
bi4 = -.0630 1 
b 32 = -.0717 , 

b 3 4 = .0195 

b 52 = -.0074 
bi 6 = -.0160 

9-88 

a n = 1*000 
a 13 = - .0705 

a 31 = -*010^- 
&33 = .036^ 
a 5 1 “ *0059 

b 02 = -.5887 
bo4 — .0720 

b 22 = .1772 

b 2 4 = .0027 

b 42 = -.0195 




TABLE III.- RITZ METHOD MODE SHAPES - Concluded 


Mode 7 

f 


a/b 

Frequency 

parameter, 

A 

Normalized mode -shape 
. coefficients 

1.0 

4.89 

: a Q1 = 0.1555 
a 03 = - • 1950 
a 21 = 1.000 
a 2 3 = .1088 

a 4 i = -.0988 

s b 01 = 0.1555 
b Q3 = -.1950 

b 2 i = 1.000 

b 23 = .1088 

b 4 i = -.0988 

1-5 

7.64 

a Ol ” - . 4^*91 
a-03 = “ * 1209 
a 2 i = 1*000 

3^3 = .04*88 

a 41 = -.0658 

b 0 i = .6498 

b 03 = -.0355 
b 21 = .5491 

b 2 3 = .1078 

t 41 = - .0590 


ro 

vo 
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Figure 1.- Method of support and excitation. 
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Figure 5.- Examples of Ritz method convergence. 







